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Abstract 

Output from the primitive equation model of Semtner and Chervin is used to examine the 
seasonal cycle in bottom pressure (pb) over the global ocean. Effects of the volume-conserving 
formulation of the model on the calculation of Pb are considered. The estimated seasonal, 
large-scale pb signals have amplitudes ranging from less than 1 cm over most of the deep ocean to 
several centimeters over shallow, boundary regions. Variability generally increases toward the 
western sides of the basins, and is also larger in some Southern Ocean regions. An oscillation 
between subtropical and higher latitudes in the North Pacific is clear. Comparison with 
barotropic simulations indicates that, on basin scales, seasonal Pb variability is related to 
barotropic dynamics and the seasonal cycle in Ekman pumping, and results from a small, net 
residual in mass divergence from the balance between Ekman and Sverdrup flows. 
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1. Introduction 


Sea level pressure has played a historical and es- 
sential role in atmospheric science, from the study 
of storm systems and numerical weather prediction 
to the characterization of climate events such as the 
El Niho/Southern Oscillation and the North Atlantic 
Oscillation. Much of this prominence has been fos- 
tered by the wide coverage and availability of baro- 
metric pressure observations. Drawing a parallel with 
the equivalent variable for the oceans, bottom pres- 
sure (pb), reveals a drastically different situation. 
There have been some focused programs that included 
pressure measurements at the deep sea floor, like 
MODE [Brown et al. , 1975] and more recently BE- 
MPEX [Luther et a/., 1990], and a few deep observa- 
tions in the Southern Ocean, mainly connected with 
studies of the Antarctic Circumpolar Current [e.g., 
Hughes and Smithson , 1996; Woodworth et al., 1996]. 
Most other available observations are from relatively 
shallow depths. Overall, though, the data coverage is 
insufficient to get even a crude estimate of the sea- 
sonal cycle over the large scale. 

In their work on the seasonal variability of the 
ocean, Gill and Niiler [1973; hereafter referred to as 
GN] assumed an ocean governed by geostrophic and 
Sverdrup dynamics in the interior, and Ekman dy- 
namics near the surface, to derive a relatively simple 
equation for p& of the form 
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where H is the ocean depth, 4> and A axe latitude and 
longitude, a and f) are the Earth’s radius and rotation 
rate, p Q is a mean ocean density, we is the Ekman 
pumping velocity, and P is the potential energy of 
the water column per unit area, defined as 


P = 



with g and p being the acceleration of gravity and 
ocean density, respectively. Variables in (1-2) repre- 
sent deviations from a time mean. The pt, equation 
should be valid for seasonal and long spatial scales 
(~1000 km or longer), away from western boundary 
regions. The two terms on the right represent, re- 
spectively, baroclinic processes associated with den- 
sity changes over the water column, and forcing by 
the wind stress curl. For given P and we fields, inte- 
gration along / jH contours (/ = 2Qsin0) from the 
eastern boundary yields pb\ the method fails in re- 


gions of closed f/H contours. GN found that, over 
the tropics, the potential energy term is important 
and baroclinic effects play more of a role, but pole- 
ward of 30° latitude, the we term dominates. Thus, 
the barotropic response to Ekman pumping at mid 
and high latitudes should play a leading role in the 
dynamics of seasonal fluctuations in p f>. Their the- 
ory remains essentially untested to date due to lack 
of appropriate data. 

The ability to characterize fluctuations in pt, and 
understand the underlying dynamics has become in- 
creasingly important in recent years. The advent of 
global altimetry has highlighted the complementary 
nature of sea level (£) and pt observations made clear 
by GN. Knowing both £ and p& could lead to more 
accurate estimates of density (or conversely steric 
height) fluctuations over the water column, and e.g. 
the possibility for differentiating between eustatic and 
steric contributions to global mean sea level trends. In 
some cases, improved estimates of surface heat fluxes 
would be possible [e.g., Wang and Koblinsky , 1997]. It 
has also become clear that the interaction of pb signals 
with topjography leads to important transfers of mo- 
mentum | between the ocean and solid Earth, on both 
short and long spatial scales [e.g. Holloway , 1992; 
Ponte a$d Rosen , 1994]. Knowledge of p& is of equal 
importance for many other current geophysical prob- 
lems, be it the determination of loading effects on the 
solid Earth [vanDam and Herring , 1994] or the inter- 
pretation of satellite geodesy measurements in terms 
of the Earth’s variable gravity field [Dong et al., 1996]. 

Given its broad geophysical relevance, there is a 
clear need for better understanding the variability in 
Pt. The Subject gains increased significance in light of 
the planned launching, in a few years, of the Gravity 
Recovery and Climate Experiment (GRACE) satel- 
lite mission, dedicated to the study of the Earth’s 
gravity f eld. Knowledge of pb field may be improved 
through indirect estimation using simple theories such 
as that of GN, or more complex modeling efforts, but 
such studies have been rare. In this paper, we at- 
tempt a characterization of the large-scale, seasonal 
cycle in pb, as simulated by the model of Semtner and 
Chervin [1992]. Some subtle issues in diagnosing pb 
from model sea level and density fields are discussed 
first. Tie seasonal cycle in pb is then analyzed and 
the unde dying dynamics are examined, in comparison 
with simulations from a simpler, barotropic model, 
and in relation to wind stress forcing. 
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2. Calculating Bottom Pressure 

In most primitive equation ocean models, pt is not 
carried as a prognostic variable and must instead be 
diagnosed from the sea level and density fields. From 
integration of the hydrostatic equation over the whole 
water column, is given by 

Pb = 9 [ pdz + p a a gp 0 C, + g [ pdz + p a (3) 

J-h J-H 

where p a is the surface atmospheric pressure. Given 
the general tendency for an inverted barometer (IB) 
sea level response to p a [Ponte, 1993; Gaspar and 
Ponte , 1997], we assume without loss of generality 
that C - C i6 + C'> i-e., IB term plus dynamic signals C' 
related to pressure, wind, or any other forcing, where 
as usual 

C b = {Pl- Pa) ( 4 ) 

gpo 

and p a is the averaged pressure over the global oceans 
(overbar will denote spatially averaged quantities through- 
out the paper). Then (3) reduces to 

— 

Pb = 9PoC +9 pdz + Pa . (5) 

J-h 

For later use, it is convenient to normalize (5) by gp Q 
to obtain 

— =C'-Cp+Ca- (6) 

gpo 

All terms have units of length, with being the steric 
height, and G being the signal in centimeters of wa- 
ter equivalent to pT Throughout the paper, bottom 
pressure will refer to the normalized pressure in (6) 
in centimeters of water. 

The diagnoses of pt, from (3) or (6) must be consis- 
tent with the model dynamics and thermodynamics, 
and the treatment of the boundary forcing. In this re- 
gard, two issues need to be taken into account. First, 
under the IB assumption, p a forcing is neglected in 
most modeling efforts. Thus, for all purposes, p a _h s 
set to zero and its effects on pt, represented by Caj 
are not considered, as they are not important dy- 
namically. For applications requiring “true^pt (as 
measured by a pressure gauge), estimates of Ca would 
have to be obtained from available atmospheric data. 

A second, more important issue in calculating pb 
has to do with the fact that most current ocean mod- 
els use a continuity equation that implies volume con- 
servation, i.e., C 7 = 0 [ Greatbatch , 1994]. However, in 
reality ocean volume (or for that matter ocean mass) 


is not strictly conserved. One way for volume changes 
to happen is through the expansion or contraction 
of water parcels associated with any density changes. 
In volume-conserving models, changes in density are 
not accompanied by corresponding volume changes. 
Thus, for example, a uniform warming of the global 
ocean, with no changes in mass involved, would still 
show up as a signal in p&, because the changes in 
are not compensated by changes in . 

To correct for the missing physics, one needs to 
consider how the ocean reacts to the unmodeled vol- 
ume changes. A simplified approach to the problem 
is discussed in detail by Greatbatch [1994] and Mel- 
lon and Ezer [1995]. Volume changes act essentially 
as a barotropic load on the ocean. Under the as- 
sumption of an equilibrium adjustment and decou- 
pling between barotropic and baroclinic processes, the 
unmodeled volume changes translate into a spatially 
homogeneous sea level signal given by 




_i 
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(7) 


(A and V are the surface area and volume of the global 
ocean), and involving the averaged density or steric 
height changes over the global ocean. Here, we shall 
take this correction to be valid for the seasonal time 
scale. At these periods, the equilibrium hypothesis 
likely holds [see Ponte , 1993]. We note, however, that 
the assumption of barotropic and baroclinic decou- 
pling needs more careful future investigation. 

Another way for effecting oceanic volume changes 
is through freshwater fluxes at the atmosphere and 
land interfaces. These effects are, however, not rep- 
resented in volume-conserving models. In fact, the 
freshwater flux condition is usually implemented as a 
restoring to observed surface salinity fields, involving 
a virtual salt flux instead of a real freshwater flux (see 
Huang [1993] for a detailed treatment of these issues). 
Again under the same assumptions used to arrive at 
(7), the missing effects of real freshwater fluxes on sea 
level can be approximated by 

c < = ~jJ A J QdtdA ' (8) 

which involves the time-integrated, spatial average of 
the freshwater flux Q (evaporation minus precipita- 
tion and runoff). Any net Q is basically spread uni- 
formly over the global oceans, with corresponding ef- 
fects on p/,. 

The sea level correction factor defined by 

(7), is easily calculated from respective model den- 
sity fields. On the other hand, calculation of ( q and 
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Co requires averaged Q and p a data, either from cli- 
matologies or atmospheric operational products, that 
are likely to have large uncertainties. Their treatment 
will not be attempted here, but there are a couple of 
points worth discussing for future reference. First, 
one should note that Cg and Ca are likely related, 
to the extent that changes in the mean atmospheric 
pressure over the oceans may involve variability in 
the atmospheric water content, through evaporation 
and precipitation processes, besides mass shifts from 
land- to ocean-covered regions and vice versa. There- 
fore, C<? and Ca effects should be considered together 
and estimated in a consistent manner. In particular, 
Ca should be calculated based on p a values that fully 
represent effects of atmospheric water content, which 
is not likely the case with at least some of the cur- 
rent operational products [van den Dool and Saha, 
1993]. Second, effects of C? and Ca on Pb m ay not be 
negligible. Operationally derived p a values point to a 
seasonal cycle in Ca with peak-to-peak amplitudes of 
0.5-1 cm [Ponte, 1993; Wunsch and Stammer, 1997]; 
similar seasonal amplitudes for C q are obtained from 
available Q climatologies. These values may seem 
small, but as discussed next, are of the same order 
of magnitude of dynamically relevant seasonal signals 
in p b . 

3. Model Estimates of the Seasonal 
Cycle in p b 

In the absence of appropriate data sets to charac- 
terize in detail the large-scale, seasonal cycle in p b , 
we examine such fields as determined from output of 
the Parallel Ocean Climate Model (POCM) of Semt - 
ner and Chervin [1992]. We use monthly averaged 
output from the so-called POCM_4B run, which is 
described in detail and analyzed in comparison with 
observations by Stammer et al [1996]. Forcing for 
this run consisted of 3-day averaged wind stresses 
calculated from European Centre for Medium-Range 
Weather Forecasts (ECMWF) twice-daily wind fields, 
ECMWF-based monthly mean heat fluxes from Bamier 
et al. [1995], and restoring conditions to Levitus et al. 
[1994] for the surface temperature and salinity fields. 
The POCM output thus represents signals in (' and 
Q p related to both mechanical wind stress forcing and 
thermodynamical forcing, but not pressure forcing. 
Given the model’s volume-conserving formulation, all 
the points discussed in the previous section are rele- 
vant in interpreting the calculated p b fields. 

The model was run initially for eight years (1987— 


1994). Monthly averaged temperature, salinity, and 
sea level fields for this period were used in (6) to calcu- 
late time series of p b at each model grid point (spacing 
of 0.4° ill longitude and 0.4° cos <j) in latitude). Given 
the focuH on large scale patterns, averages over 4°x4° 
boxes were constructed and used in the analyses. In 
addition, to minimize any effects of transient adjust- 
ment of the model fields to the new boundary forcing 
over the initial period, the analysis is restricted to the 
last six years of output (1989-1994). 

Beforb analyzing the spatial character of the sea- 
sonal variability, the influence of volume effects as- 
sociated with 6£p are examined. Figure 1 shows the 
time series of 6£ py calculated by integrating the den- 
sity field! over the volume of the model as in (7), with 
p 0 = 1.0$ g/cm 3 and A = 3.47 x 10 18 cm 2 , the surface 
area of tihe model ocean. There is a long term trend 
of -2 cih over the six year period (approximately 3 
mm/yeat), which is fairly linear over the first four 
years and tends to flatten out in the last two years. 
Models are known to have a tendency to drift with 
rime, due to slow diffusive adjustment, incompatibili- 
ties in boundary forcing and model physics, etc. (see 
Stammer et al. [1996] for a brief discussion of drifts in 
POCM). and the implicit decrease in global mean sea 
level in Fig. 1 (model ocean is getting denser) seems 
too large to be realistic. Examination of these trends 
should be an important consideration for studies of 
interannual and lower frequency variability in p b , but 
need not concern us in the study of the seasonal cycle. 

Superposed on the long term trend in Fig. 1, there 
is a clear seasonal cycle in 6Q P with peak-to-peak am- 
plitudes of approximately 0.3-0. 4 cm, and minimum 
occurring around August. This signal is a measure 
of hemispheric asymetries in the ocean-land distri- 
bution aind the seasonal cycle in steric height. Mel - 
lor and Ezer [1995] have calculated 8Q P from avail- 
able hydrographic data and estimate a seasonal cy- 
cle with a peak-to-peak amplitude of approximately 
0.8-1 cm, and minimum in July (see their Fig. 11). 
The tendency of POCM to underestimate the sea- 
sonal cycle in steric height, apparently due to lack 
of proper heat uptake by the mixed layer, results in 
a weaker seasonal signal in 8( p , and the phase shift 
by approximately one month with respect to the ob- 
servatior s may be partly due to phase errors of the 
forcing ECMWF heat fluxes [Stammer et al., 1996]. 
The diffe rent 6£ p may of course be also due to the cli- 
matological, sparse hydrographic records used in Mel- 
lor and Ezer’s calculation, as opposed to the 6-year, 
dense sampling of the model. 
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The pf, fields analyzed here include the model cal- 
culated term shown in Fig. 1. After removing 
means and trends from each pb time series, values 
were averaged for each season (December-February, 
March-May, June- August, and September-November). 
The corresponding mean seasonal variability in pb 
over the 6-year period of model output is shown in 
Fig. 2. Most of the deep ocean, including tropical lat- 
itudes, exhibit pb fluctuations on the order of 1 cm 
or smaller. (Note that these amplitudes are of the 
same order of magnitude of the seasonal cycle in 
the effects of 5( p should thus be taken into account 
in the study of pb , in contrast with their less impor- 
tant role for the variability in £■) Larger variability 
is seen over the Southern Ocean, the western North 
Pacific, and near boundary regions, particularly over 
shallow depths (e.g., North Sea, Bering Sea, Sea of 
Okhotsk, South China Sea), where amplitudes of a 
few centimeters are attained. 

The large, basin scale character of seasonal pb sig- 
nals is striking. The Atlantic, Indian, and to some 
extent the Southern Ocean tend to show single sign 
patterns. In northern winter, most of the Atlantic 
and Southern Oceans show a positive mass anomaly, 
with the Pacific and Indian Oceans having a negative 
anomaly. The pattern is reversed in northern sum- 
mer. Besides the apparent seasonal interbasin mass 
exchange, there is a clear oscillation in the North Pa- 
cific, with a node at approximately 30°N, which sug- 
gests seasonal strengthening of the depth-averaged 
gyre circulation in northern winter and spring, and 
mass exchanges between subtropical and higher lati- 
tudes. 

To assess the seasonal variability in Fig. 2 as a 
function of frequency, Fig. 3 shows the power spec- 
trum of pt, calculated by averaging the periodogram 
over all grid boxes in the model. Each time series was 
detrended and tapered using a cosine bell function 
over 10% of the points on each end, prior to Fourier 
transformation. Detrending does decrease the power 
slightly at the first two harmonics (6 and 3 years), but 
has no effect on higher harmonics. The pb spectrum 
shows a large dominant peak at the annual period, 
and a weaker peak at the semiannual period. Besides 
these peaks, the spectrum exhibits a dependence on 
frequency close to u; -0,5 . The important contribution 
of the annual time scale to the seasonal variability in 
Pb , seen in Fig. 2, is clear. 

Given the focus of many studies on the annual pe- 
riod, we show in Plate 1 the amplitude and phase of 
the annual cycle in pb, calculated by simply Fourier 


transforming each time series after detrending. Most 
of the features discussed in relation to Fig. 2 are 
present in Plate 1. The annual cycle has amplitudes 
smaller than 1 cm over most of the global oceans, 
with higher values (1-3 cm) over the Southern Ocean 
and the western North Pacific, and over most of the 
shallow coastal regions in the model. The phase ex- 
hibits, for the most part, a large scale pattern: the 
lack of a distinct phase opposition between the North- 
ern and Southern Hemispheres, observed in the case 
of C and Cp, is most striking. Maximum values in 
the Atlantic and Pacific Oceans occur mostly in the 
northern winter (Dec-Feb) and summer (Jun-Aug) 
seasons, respectively. Thus, the Atlantic and the Pa- 
cific appear to fluctuate approximately out of phase, 
as noted by Stammer et ai [1996], based on the av- 
eraged pb field over the Atlantic and Pacific Oceans 
(see their Fig. 15). 

Before attempting to interpret the seasonal pb sig- 
nals, it is also useful to examine the relation between 
Pb and by calculating 


(gpo) 

no 


( 9 ) 


which is essentially the admittance between pb and 
with T denoting Fourier transform. In physical 
terms, A quantifies how much of a bottom pressure 
signal results from a sea level load at the surface. 
From (6), one expects |.4| < 1 if Q p and £' are cor- 
related. Likewise, phases near zero indicate a sim- 
ple in-phase relationship among and Cp* The 

amplitude and phase of A for the annual period are 
shown in Plate 2. Variability in pb is weak relative 
to with |*4| < 0.4 over most of the oceans, away 
from shallow regions and the Southern Ocean. The 
small pb amplitudes are expected as a result of the 
well-known correlation between and Cp- There is, 
however, extensive regions, particularly in the South- 
ern Ocean where \A\ > 0.4. Over most shallow coastal 
regions, strong correspondence between pb and £' is 
indicated by amplitudes near unity and phases near 
zero. These results most likely indicate the essen- 
tially barotropic nature of variability in these regions. 
A tendency for nearly in-phase behavior is clear in 
the Southern Ocean and the North Pacific (north of 
~ 30°), but phases span the whole range. In par- 
ticular, the approximately out-of-phase behavior over 
equatorial Pacific regions is worth noting, as indica- 
tive of the presence of first (or more precisely odd) 
baroclinic mode variability [e.g., Brown et al., 1975]. 
The regional dependence of the annual admittance in 
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Plate 2 suggests a complex relation between p b and 
C and the difficulty in extracting information about 
pt from the sea level fields alone. 

4. Dynamics and Forcing of Seasonal 
Variability in p b 

As discussed in the Introduction, the importance 
of barotropic dynamics and wind stress curl forcing 
for large-scale, seasonal variability in p b at mid and 
high latitudes was pointed out by GN. Model results 
in Fig. 2 are qualitatively similar to the North At- 
lantic and North Pacific estimates of GN (see their 
Fig. 3), but more general evidence for the impor- 
tance of barotropic dynamics to the seasonal p b vari- 
ability can be obtained by comparing the POCM re- 
sults to those from a recent run of the near-global, 
shallow-water model of Ponte [1993]. For details 
about model formulation, resolution, friction, etc., 
the reader should consult Ponte [1993] and Caspar 
and Ponte [1997]. The available model output is from 
a run originally performed for other purposes and 
forced by four-times daily reanalysis wind stress and 
pressure fields from the National Centers for Environ- 
mental Prediction (NCEP) for the years 1993-1995, 
with a spinup period taken to be the last two months 
in 1992. Thus, C represents the effects of both pres- 
sure and wind stress, but the latter dominates at the 
seasonal time scale [Ponte, 1994]. For the barotropic 
model, there is of course no role for thermodynam- 
ics, Cp is assumed constant, and calculation of p b is 
equivalent to that of 

From the four- times daily model output, seasonal 
means of (' for the months Dec-Feb, Mar-May, J un- 
Aug, and Sept-Nov were calculated. The resulting 
seasonal variability, based on years 1993-94, is shown 
in Fig. 4 for comparison with POCM results derived 
for the same 1993-94 period and shown in Fig. 5. 
The comparison is qualitative, as different operational 
products were used for the forcing fields in the two 
models, in addition to differences in resolution, to- 
pography, coastal geometry, etc. Nevertheless, be- 
sides a noticeable tendency for stronger signals in the 
barotropic model, the similarities between the POCM 
and barotropic results are clear. Correspondence be- 
tween patterns of positive and negative p b anomalies, 
and regions of larger and smaller amplitudes, is qual- 
itatively good, especially for the winter and summer 
seasons. The western intensified, North Pacific os- 
cillation between subtropical and higher latitudes is 
particularly well repraduced in the barotropic results. 


Some of the characteristics of the p b seasonal vari- 
ability in Figs. 4 and 5 can be traced to the na- 
ture of the wind stress curl and the corresponding 
Ekman pumping fields in (1). Figure 6 shows the sea- 
sonal variability of the forcing we term in (1), calcu- 
lated from monthly averages of the NCEP wind stress 
fields for the years 1993-94 that were used to force the 
barotropic model. Derivatives of the wind stress in- 
volved in the calculation were evaluated at every grid 
point using centered differences; the calculation was 
confined to latitudes poleward of 5°, as Ekman the- 
ory formally breaks down at the equator. Without 
attemptihg a full integration along f /H contours as 
done by GN, we note that the seasonal forcing sig- 
nals are generally stronger and more zonally homoge- 
neous in the North Pacific than in the North Atlantic. 
These two factors provide for the stronger p b variabil- 
ity and the clear dipole pattern in the North Pacific, 
as compared to the North Atlantic (the longer extent 
of the Pacific may also lead to stronger variability to- 
ward the western boundary). More generally, there is 
qualitative correspondence between the positive and 
negative large scale patterns of the we and p b fields, 
respectively, and vice versa, which is consistent with 
(1). The| correspondence is particularly good for sea- 
sons and basins for which the we fields are more spa- 
tially homogeneous. 

The relation between p b and we is also apparent 
when examining their frequency spectra. It can be 
easily shown from (1) that, if the spectrum of we is 
separable in frequency and wavenumber, i.e., 

j 4>(K,u) ~ (t>i(K)(f>2{uj) ( 10 ) 

then the frequency dependency of the p b spectrum is 
simply given by Given that the spectrum of 

we is essentially that of the wind curl, the approx- 
imation in (10) is a reasonable one at low frequen- 
cies [e.g., Willebrand et al. , 1980]. We have calcu- 
lated the periodogram of the we term in (1) at every 
grid point (again excluding latitudes within 5° of the 
equator) and averaged them together to obtain the 
frequency spectrum shown in Fig. 7. The frequency 
dependency of that spectrum is remarkably similar to 
the sped rum ofp*, shown in Fig. 3, both in its general 
backgroi nd shape (close to u“° 5 ) and in the relative 
height of the peaks at annual and semiannual periods. 

Seasonal fluctuations in p b field, highlighted in 
Fig. 2 and Plate 1, represent large-scale divergence 
or convergence of mass brought about by the season- 
ally var^ ing flow field. Horizontal oceanic motions 
are largely nondivergent, however, and only a small 
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divergence is involved in effecting pb changes. For ex- 
ample, to accumulate 1 cm of water over a 10 4 xl0 4 
kilometer area in six months, one needs only a net 
transport of ~6xl0 4 m 3 s -1 , much smaller than typi- 
cal horizontal transports in the ocean. To explore the 
dynamics behind the divergent mass fluxes leading to 
the pb signals, it is useful to consider a more general 
form of (1), which accounts for more terms in the di- 
vergence field. The derivation is given in Appendix A 
using Cartesian coordinates and neglecting nonlinear 
and dissipation effects for simplicity. The equation 
for pb is 





— p 0 Q + O 


lj 

7 


terms. 


(ii) 


As pointed out in Appendix A, (11) associates local 
mass divergence, represented by pb t , with four dif- 
ferent processes: geostrophic flow divergence due to 
variable / and H, which includes the baroclinic term 
P x ; Ekman flow divergence forced by the wind stress; 
mass fluxes at the surface due to Q\ and ageostrophic 
effects of 0[tj/f] relative to the geostrophic term. 

Introducing a spatial scale of motion L, pb t term is 
0[(lj/ f)(L/L,R) 2 ) relative to geostrophic term, where 
Lr — y/gU / f is the external Rossby radius (~ 2000 
km). Seasonal Q forcing is generally negligible com- 
pared to the Ekman terms. For seasonal time scale 
and away from the equator (u )/ f <C 1), and for typ- 
ical L ~ 1000 km, one has the traditional quasi- 
steady Sverdrup balance, as in (1), with wind-driven 
Ekman divergence being for the most part balanced 
by geostrophic flow divergence. The net residual di- 
vergence is associated with small geostrophic flow 
correction to the Sverdrup solution, at scales larger 
than Lr . Basin-scale signals discussed in Fig. 2 are 
likely dominated by these processes. At scales shorter 
than Lr, divergence from 0[u/f) ageostrophic effects 
should also become important, as well as dissipation 
and nonlinear effects not treated in (11). 


5, Final remarks 

It is difficult to assess the realism of the model 
p b fields in Fig. 2 in the absence of global pressure 
records. From the available point measurements over 
the deep ocean, only a few are longer than one year 
[Woodworth et ai, 1996; Mofjeld et ai , 1996]. The 
observations just to the east of the Drake Passage, an- 
alyzed by Woodworth et al. [1996], show fluctuations 


over a variety of time scales; amplitudes of the sea- 
sonal cycle are not provided, but visual inspection of 
their Figs. 5a, 6, 7, 8a is not inconsistent with ampli- 
tudes at the centimeter level. The records of Mofjeld 
et ai [1996], taken near 46°N, 130°W in the North 
Pacific, show substantial contamination from instru- 
ment drift after 30-day running means are applied to 
them (Mofjeld, personnal communication, 1997), but 
all indications point to seasonal amplitudes with or- 
der of magnitude similar to the records from Wood- 
worth et al. Not much can be said quantitatively 
about the agreement or disagreement between model 
pb and these records, except to note that observed val- 
ues seem larger than (but of the same order of mag- 
nitude of) the POCM results. Point measurements 
may of course relate to local processes not relevant 
to large-scale variability examined in the model. Di- 
rect, more quantitative evaluation of model estimates 
of pb against observations must await the launching 
of the gravity missions currently being planned. Con- 
versely, the seasonal fluctuations in the mass field over 
the oceans estimated here provide a crude measure of 
the oceanic signals that will need to be resolved by 
the planned gravity missions. 

Besides comparisons with point measurements, there 
are other more indirect methods by which the model 
Pb fields can be evaluated. Given the importance of 
Pb variability in the excitation of Earth rotation sig- 
nals, particularly seasonal polar motion [Ponte et al . , 
1998], and also its effects on the variable gravity field, 
as observed by satellite tracking methods [Dong et 
ai, 1996], estimation of these oceanic effects based 
on the model pb fields and comparison with other rel- 
evant geophysical observations on Earth orientation 
and gravity field provide an independent consistency 
check on the pb fields, as shown in the case of a dif- 
ferent ocean model by Ponte et ai [1998]. Analysis 
of the pb fields from the POCM model in this broader 
geophysical context is currently being pursued by this 
and other investigators. 

In all these comparisons, one must keep in mind 
that signals related to and may contribute to 
the local pb measurements, and may also have a sig- 
nificant projection onto the integrated measurements 
represented by Earth rotation and large-scale grav- 
ity parameters. Given the crude estimates of ( a and 
Q available, and depending on their relative phase, 
these effects may impact the estimates of the seasonal 
cycle in regions with pb signals at millimeter ampli- 
tudes. Ideally, one would like to have pressure forcing 
and real freshwater fluxes implemented in circulation 
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models, which would permit verification of the equi- 
librium response assumption at seasonal time scales 
and examination of higher frequency variability where 
that assumption becomes less useful [Ponte , 1993]. 

Finally, regarding the inferrence of surface heat 
fluxes or heat content in the water column from al- 
timetry sea level records, as attempted recently by 
Wang and Koblinsky [1997], Chambers et al [1997], 
and also White and Tai [1995], Plate 2 provides a 
crude estimate of the error induced by neglecting con- 
tributions to C which are truly related to pt> and not 
for the annual cycle. Depending on the region, es- 
timates may easily be in error by 20 to 40%. However 
large these values may seem, the sea level-based esti- 
mates are still meaningful given the large uncertain- 
ties in the surface heat fluxes and the lack of sufficient 
hydrographic data sampling to estimate heat content. 
The possibility of improving those estimates by using 
model-derived pb fields is, nevertheless, noted. 


Appendix A: Derivation of General 
Equation for p ^ 

Starting from the linearized horizontal momentum 
and cont inuity equations, vertical integration over the 
full depth H leads to 

/ 1 f° , T x 

v ' +ln ‘-JJiL l , p ’ iz + W* (M) 

Ct + ( Hu) x + ( Hv)y = C p,-Q 

Here, u and v are depth-averaged flows, £ is sea level, 
C p is stdric height, and Q denotes freshwater flux 
(evaporation minus precipitation and runoff). We 
work with Cartesian coordinates and neglect bot- 
tom and j lateral dissipation effects for simplicity. The 
ocean is acted upon by surface wind stresses r x,y . 

To write (Al) in terms of pi,, one needs to work on 
the pressure integrals. Moving the derivative outside 
the integral, and using the fact that 


p(z) = Pb - 9 
the integrals become 


f 

J-H 


pdz 


f H P x , v i 2 = ^HPb - g j dz J pdz’^j -p b H x<v 

The double integral can be integrated by parts and 
after carrying out the derivatives, we have simply 


/ 


Px.ydz = Hp biy + P, 


H 


*.V 


(A 2 ) 


with P as defined in (2). 

With -he use of (A2) and the hydrostatic relation 


— = c 

gpo 


C p 


(A3) 


(effects cf p a are neglected here), we can write (Al) 
as 

f Pb x Px , T x 

u t — fv = 

Po Hpo Hp 0 


Vt -b fu — — 


^ ^ — (A4) 

p 0 Hp 0 H Po ^ ; 


Pb< 

9Po 


+ (Hu) x + (Hv) y = -Q 
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These equations are very similar to the linear shallow- 
water equations forced by wind stress and freshwater 
fluxes, with replacing £ and extra forcing terms 
P x ,y associated with baroclinic effects. With time de- 
pendence of the form e luJt , it is straightforward to 
solve the momentum equations for u,v and substi- 
tute in the continuity equation to get a very general, 
but complicated, equation for For seasonal time 
scales and away from the equator, i.e., tu /, and 
without explicitly treating 0[u / /] terms arising from 
Ut,Vt terms in (A4), one arrives at a simpler equation: 

Pb t r ( fy p 

9 \ Pb, f) P X ~ 

(jj'j ~ (j^) ~ Po ® + ° J terms (^ 5 ) 

If one partitions the flow Held in the usual geostrophic 
and Ekman components in (A4), it is easy to see that 
the wind stress terms in (A5) represent the divergence 
of the Ekman transports, whereas the Jacobian and P 
terms represent the divergence of geostrophic trans- 
ports. Ageostrophic effects due to terms in mo- 
mentum equation are represented in 0[u>/ f] terms. 
Any imbalance of mass fluxes due to these Ekman, 
geostrophic, and ageostrophic flows, and to surface 
mass inputs ( Q term), lead to variability in 
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Figure 1. Time series of £ p (cm) for the period January 1989-December 1994 calculated from POCM output as 
described in the text. Origin of the time series is arbitrarily set to zero. 

Figure 2. Anomalies in p b averaged over the months of December-February, March-May, June- August, and 
September-November. Contour interval is 0.5 cm; negative values are lightly shaded. 

Figure 3. Power spectrum obtained by averaging individual periodograins of p b (no area weighting) calculated at 
all model points. No averaging in frequency is applied. 

Figure 4. Anomalies in p b as in Fig. 2, but calculated from the barotropic model based on two years of output 
(1993-1994). Countour intervals are 1 cm; negative values are lightly shaded. Dark shading represents land points 
in the model. 

Figure 5. Anomalies in p b as in Fig. 2, but calculated from the POCM model based on only two years (1993-1994). 
Countour intervals are 1 cm; negative values are lightly shaded. 

Figure 6. Seasonal anomalies in the Ekman pumping forcing term in (1), based on wind values for 1993-94. 
Values within 5° of the equator were not calculated. Countour interval is 1(T 4 g s -2 ; negative values are lightly 
shaded. The land mask is that of the barotropic model in Fig. 4. 

Figure 7. Power spectrum of the Ekman pumping term in (1), calculated by averaging periodograms at all points 
for which data is shown in Fig. 6. No averaging in frequency is applied. 


Plate 1. Amplitude (in cm) and phase of the annual cycle in p 6 . Note the changing contour intervals for the 
amplitude. Phases are plotted every other point and indicated by the direction of the line segments (origin denoted 
by the dots); 0° corresponds to eastward direction, and phase increases anticlockwise. Values between -45° and 
45° correspond approximately to a maximum in northern winter season (Dec Feb), between 45° and 135° to a 
maximum in northern spring, and so on. 


Plate 2. Amplitude and phase of the admittance between p b and (' for the annual harmonic. Phase is plotted as 
in Plate 1. Line segments pointing eastward (westward) imply in-phase (out-of-phase) relation between p b and 
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